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ABSTRACT 

We present measurements of the column densities of interstellar D I, O I, N I, and H2 made with 
the Far Ultraviolet Spectroscopic Explorer (FUSE), and of H I made with the International Ultraviolet 
Explorer (IUE) toward the sdO star LSE 44 [(l,b) = (313?37, +13? 49); d = 554 ± 66 pc; z = +129 ± 15 
pc] . This target is among the seven most distant Galactic sight lines for which these abundance ratios 
have been measured. The H I column density was estimated by fitting the damping wings of interstellar 
Lya. The column densities of the remaining species were determined with profile fitting analyses, and 
supplemented with curve of growth analyses for O I and H 2 . We find log N(D I) = 15.87 ± 0.08, log 
N(0 I) = 17.57±o i5, log 2V(N I) = 16.43 ± 0.14, and log N(R I) = 20.52±°j° (all errors 2a). This 
implies D/H = (2.24±[]H) x 10" 5 , D/O = (1.991^?) X 10~ 2 , D/N = (2.75±S;||) x 10" 1 , and O/H = 
(1.131 

an) x 10 3 - Of the most distant Galactic sight lines for which the deuterium abundance has been 
measured LSE 44 is one of the few with D/H higher than the Local Bubble value, but D/O toward all 
these targets is below the Local Bubble value and more uniform than the D/H distribution. 

Subject headings: Cosmology: Observations - ISM: abundances - Ultraviolet: ISM - stars: Individual 
(LSE 44) 



1. INTRODUCTION 

Precise measurements of primordial abundances of the 
light elements deuterium (D), 3 He, 4 He, and 7 Li relative 
to hydrogen have been a goal of astronomers for many 
years. In the standard Big Bang nucleosynthesis model, 
these quantities are related in a straightforward way to 
the baryon-to-photon ratio in the early universe, from 
which f2e, the fraction of the critical density contributed 
by baryons, may be determined (Boesgard & Steigman 
1985). Since deuterium is easily destroyed in stellar interi- 
ors (astration), and no significant production mechanisms 
have been identified (Reeves et al. 1973; Epstein, Lat- 
timer, & Schramm 1976), the D/H ratio is expected to 
monotonically decrease with time. 

Precise measurements of D/H in the local interstellar 
medium (ISM) were first made using Copernicus (Roger- 
son & York 1973), HST (Linsky et al. 1995), and IMAPS 
(Jenkins et al. 1999; Sonneborn et al. 2000). In the last 
several years many additional measurements of D /H, D/O, 
D/N, and related species, have been made using FUSE 
(Moos et al. 2002; Hebrard & Moos 2003; Wood et al. 
2004). These local ISM measurements represent a lower 
limit to the primordial value. When corrected for the ef- 
fects of astration (Tosi et al. 1998), the results should 
be comparable to those obtained in low-metallicity, high 
redshift Lya clouds (Crighton et al. 2004; Kirkman et al. 
2003), which should be very nearly the primordial value 
itself. 

As data on more sight lines has accumulated the situ- 
ation is not as simple as originally envisioned. At high 
redshift, measurements of D/H vary from (1.6 — 4) x 10 -5 
(Crighton et al. 2004; Pettini & Bowen 2001), despite 



the fact that the metallicity in all these environments is 
low enough to suggest that D/H should match the pri- 
mordial value. More locally, it is now accepted that for 
sight lines within the Local Bubble, a region of hot, low- 
density gas (Sfeir et al. 1999) within a distance of 100—150 
pc and log JV(H I)< 19.2, in which the Sun lies, D/H 
is nearly constant, but the value is uncertain. Based 
on measurements toward 16 targets, Wood et al. (2004) 
state that (D/H) lb = 1-56 ± 0.04 x 10 -5 . Alternatively, 
(D/H) lb may be inferred from measurements of D/O and 
O/H. Hebrard & Moos (2003) discuss the advantages of 
this approach, including potential reduction of systematic 
errors associated with measuring N(D I)and N(H I) di- 
rectly, which differ by five orders of magnitude. Their 
approach relies on the observed high spatial uniformity of 
O/H (Meyer, Jura, & Cardelli 1998; Meyer 2001; Andre 
et al. 2003). The result is (D/H) LS = 1.32 + 0.08 x 10~ 5 , 
about 2a different than the Wood et al. value. 

At distances greater than ~ 500 pc, or more precisely, 
at log N(H I);> 20.5, the few sight lines measured have 
lower D/H. Based on a suggestion originally made by Jura 
(1982) and extended by Draine (2004a, 2004b) that un- 
der the right thermodynamic conditions deuterium can be 
depleted onto dust grains, Wood et al. (2004) have in- 
terpreted D/H variations as being due to various levels 
of depletion and, over the longer sight lines, the effects 
of averaging the D/H concentration over multiple clouds 
along the sight line. In this scenario, the low value at 
large distances represents the gas phase abundance in the 
local disk, (D/R) LDg = (0.85 + 0.09) x 10" 5 , and one must 
include the depleted deuterium to determine the total deu- 
terium abundance. An alternate explanation is offered by 
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Hebrard & Moos (2003). Based on D/O, D/N, and D/H 
measurements, they suggest that the low gas-phase value 
observed at high column densities and distances is the true 
Galactic disk value, and there is not a significant reservoir 
of deuterium trapped in dust grains. D/H has an unusu- 
ally high value in the Local Bubble if this interpretation is 
correct. The reason for this high value is not understood. 
One possibility is the infall of primordial, deuterium-rich 
material, and another is non-primordial deuterium pro- 
duction, such as in stellar flares. This is discussed further 
in §5. 

To test this hypothesis it is necessary to make additional 
measurements of D I, H I, O I, N I, H2, and additional 
species that might yield information about depletion and 
physical conditions, toward targets with log N(H 1)^ 20.5 
In this paper we present the results of an analysis of the 
sight line toward LSE 44, an sdO star located at a distance 
of 554 ± 66 pc and well within this column density regime. 
We derive D I, O I, N I, and H2 column densities with data 
obtained from the Far Ultraviolet Spectroscopic Explorer 
(FUSE) and estimate the H I column density with data 
from the International Ultraviolet Explorer (IUE). 

In §2 we present a summary of the observations and a 
description of the data reduction processes. In §3 the stel- 
lar properties of LSE 44 are discussed. The details of the 
analysis of the D I, O I, N I, and H2 column densities are 
given in §4, and the H I analysis is presented in §5. In §6 
we discuss the results of this study. 

2. OBSERVATIONS AND DATA PROCESSING 

The FUSE instrument consists of four co-aligned spec- 
trograph channels designated LiFl, LiF2, SiCl, and SiC2, 
named for the coatings on their optics, which were selected 
to optimize throughput. Each channel illuminates a pair 
of microchannel plate segments, labeled A and B. The LiF 
channels have a bandpass of approximately 1000-1187 A 
and the SiC channels approximately 910-1090 A. Thus, 
each observed spectral line may appear in the spectra from 
multiple channels. The FUSE resolution is approximately 
20 km s _1 (FWHM), but it varies slightly between chan- 
nels, and within a channel, as a function of wavelength. 
A detailed description of the FUSE instrumentation and 
performance is given by Sahnow et al. (2000). 

LSE 44 was observed with FUSE for a total of 86 ksec 
under programs P2051602 - P2051609. A total of 119 ex- 
posures are in the MAST archive, obtained with various 
levels of alignment success for the separate channels, and 
some were obtained when one or both detectors were not 
at full voltage. The observation log is shown in Table 1, 
and the total exposure time by channel is shown in Table 2. 
All data were obtained in time-tag mode with focal-plane 
offsets (FP-SPLITS) in the MDRS (20" x 4") aperture. 

The data were reduced using CALFUSE pipeline ver- 
sion 3.0.6. We used the default pixel size of 0.013A. 
The signal-to-noise ratio per pixel in the continuum re- 
gion around 928A for both SiClB and SiC2A spectra is 
typically about 3.5 in individual, well-aligned exposures of 
~ 900s, the maximum duration of any exposure. The in- 
dividual 1-dimensional spectra were co-added to form the 
final spectrum for each channel and detector segment sepa- 
rately after removing the relative shifts between individual 
spectra on the detector caused by image and grating mo- 
tion (Sahnow et al. 2000). The shifts were determined 



by cross-correlating the individual spectra over a limited 
wavelength range which contained prominent spectral fea- 
tures but no airglow lines. Typical shifts were <; 4 pix- 
els, or < 0.050A. The S/N ratios in the summed spectra 
are approximately 26 and 27, respectively, for the SiClB 
and SiC2A channels. Narrow interstellar absorption fea- 
tures in the summed spectra do not exhibit unexpected 
broadening, indicating that the summing process did not 
smear the data. Although FUSE observations of H I and 
O I lines are sometimes contaminated by airglow emission, 
data obtained during spacecraft orbital night showed that 
this was unimportant in our analysis, so the combined day 
and night data were used. 

FigureQ]shows the co-added SiClB spectrum of LSE 44, 
which covers the wavelength range 915 - 990A. D I, O I, 
N I, and H2 lines arising in the ISM are identified. 

The H I column density along this line of sight was es- 
timated using the Lya profile from the single high resolu- 
tion spectrum of LSE 44 taken with IUE. The 16.5 ksec 
observation was made on 29 August 1986 (see Table 1), 
and reduced using IUESIPS. The analysis of this data is 
discussed in more detail in §4.4. 

3. STELLAR PROPERTIES OF LSE 44 

In his survey of subluminous O stars, Drilling (1983) 
obtained photometric and spectroscopic data for LSE 44. 
His UBV photometry gives V = 12.45, B -V = -0.24, 
and U — B = —1.18 (Table 3), and indicates that the 
star is relatively hot. He describes the optical spectrum 
of LSE 44 by saying that the only lines present are the 
He II A4686 line and the Balmer series. He adds that 
the appearance of these lines is very similar to those in 
the spectrum of the sdOB star Feige 110, but the Balmer 
lines are stronger and the He II A4686 line is weaker in 
the spectrum of LSE 44 (see, e.g., Heber 1984; Friedman 
et al. 2002). In order to measure the atmospheric pa- 
rameters of LSE 44, D. Kilkenny (2004, private commu- 
nication) obtained an optical spectrum of LSE 44 at the 
South African Astronomical Observatory with the 1.9-m 
telescope and CCD spectrograph. The optical spectrum 
covers the wavelength range 3400 - 5400 A with a resolu- 
tion of about ~ 3 A. The normalized spectrum is plotted 
in Figure |21 As reported by Drilling (1983), our optical 
spectrum shows broad Balmer lines starting from H/3 to 
H6* and the He II A4686 line, which is the only helium line 
detected in this spectrum. The Ca II K A3935 line is also 
detected in the optical spectrum. 

The atmospheric parameters of LSE 44 were obtained 
by fitting the optical spectrum illustrated in Figurc[2]with 
a grid of synthetic spectra and by using the Marquardt 
method (see, Bevington & Robinson 1992). This method 
optimizes the effective temperature, gravity, and helium 
abundance in order to obtain a model that best matches 
the observed spectrum. We computed a grid of NLTE 
H/He atmosphere models with 20,000 K < T cS < 80,000 K 
in steps of 2,000 K, for different values of the surface grav- 
ity 4.8 < log<? < 6.4 in steps of 0.2 dex, and for differ- 
ent helium abundances -4.0 < log N(Re)/N(H) < 0.0 
in steps of 0.5 dex. The models were computed with 
the stellar atmosphere codes TLUSTY and SYNSPEC 
(see, e.g., Hubeny & Lanz 1995). Figure |21 shows our 
best fit for the optical spectrum of LSE 44. We ob- 
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tain T eff = 38,700±1,000 K, log 5 = 5.5 ± 0.1, and 
logiV(He)/7V(H) = -2.8±0.1. We used these atmospheric 
parameters to compute the stellar Lya line profile in order 
to measure the H I column density toward LSE 44. (see 
§4.4). 

The evolutionary status and physical properties of 
LSE 44 can be estimated by comparing its position in 
a T e f{-g diagram with evolutionary sequences. Figure [3] 
shows the position of LSE 44 in such a diagram with 
the positions of a sample of hot subdwarf B stars (sdB). 
It also shows post-extreme horizontal branch (post-EHB) 
evolutionary sequences that were computed by B. Dorman 
(1999, private communication). The position of LSE 44 in 
the Teff-5 diagram suggests that it is an sdB star that has 
evolved from the EHB. Even though LSE 44 is classified 
as an sdO star based on its optical spectrum, it relates to 
the sdB stars. sdB stars have effective temperatures and 
gravities in the ranges 24,000 K <; T cff <; 40,000 K and 
5.0 <J logg <J 6.2, and helium abundances that are typi- 
cally a factor of ten smaller than solar (see, e.g., Saffer et 
al. 1994). The evolutionary sequences in Figure |3| repre- 
sent the evolution of typical sdB stars that have helium- 
burning cores of ~ 0.477 M Q and small hydrogen envelopes 
-Wenv <^ 0.0045 M Q . Each evolutionary track corresponds 
to models with about the same helium core mass but with 
different hydrogen envelope masses. In these models the 
hydrogen envelopes are too small to ignite the hydrogen 
burning shell, and therefore the stars appear subluminous, 
hot, and compact. In Figure [3] the hatched region corre- 
sponds to the zero-age EHB (ZAEHB). In this way, an 
EHB star evolves from the ZAEHB and moves first to the 
right and then to the left toward the white dwarf region 
without reaching the asymptotic giant branch (AGB). A 
star spends ~ 10 s yr on the EHB and ~ 10 7 yr off the 
EHB en route to the white dwarf region. 

The spectral type sdO is a spectroscopic class that in- 
cludes stars showing strong Balmer lines and He II lines. 
Because of its high effective temperature and low helium 
abundance, LSE 44 appears as an sdO star. The evolu- 
tionary tracks surrounding the position of LSE 44 show 
that stars in this portion of the T e g-g diagram have left 
the helium core burning phase and undergo a helium shell 
burning phase. One can assume that the mass of LSE 44 
is ~ 0.4775 M and estimate its distance by using the 
expression 
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where G is the gravitational constant, M is the mass of the 
star, g is the gravitational acceleration, H v is the Edding- 
ton flux weighted by the Johnson passband V of Bessel 
(1990), Fo t v is the average absolute flux of Vega at V as 
specified in Heber et al. (1984), V is the apparent visual 
magnitude, and Ay is the interstellar absorption. By using 
the best model fit parameters and the passbands of Bessel 
(1990) we computed the intrinsic color index (B — V)o and 
obtained the color excess E(B — V) = 0.05 ± 0.03, which 
gives A v = 0.155 ± 0.093 when using R v = 3.1. The 
distance to LSE 44 is then 554 ± 66 pc. 

The FUSE spectrum of LSE 44 shows many photo- 
spheric and interstellar absorption lines. The strongest 
photospheric lines are the Lyman series that start at Ly/3 



and go up to lines close to the series limit. Because of the 
high gravity the photospheric Lyman lines merge together 
before reaching the series limit, and the remaining Lyman 
lines close to this limit are mainly due to interstellar H I. 
The other strong photospheric lines are He II A1084, C III 
AA1175, N III AA1184, N IV AA923 and A955, Si IV AA1126 
and A1066, P V AA1123, S IV AA1070, and S VI AA937. 
There are also many fainter unidentified stellar lines that 
depress the continuum. For instance, Figure 0] shows a 
comparison between the FUSE spectra of LSE 44 and the 
sdOB star Feige 110 (T eg = 42,300 K and \ogg = 5.95; 
Friedman et al. 2002). The figure illustrates that both 
stars have many lines in common, but because of the lack 
of atomic data, we cannot identify these absorption lines. 
The Si IV lines at 1122.49 A and 1128.33 A and the P V 
line at 1128.008 A are practically the only stellar lines that 
can be identified in these portions of the FUSE spectrum. 
This example illustrates clearly that blends of unidentified 
stellar lines with interstellar lines could lead to systematic 
errors when measuring interstellar column densities. 

Table 3 summarizes the important properties of LSE 44 
and the sight line to this star. 

4. INTERSTELLAR COLUMN DENSITY ANALYSES 

Two techniques have been used to estimate the column 
densities in this analysis. The first involves directly mea- 
suring the equivalent widths of absorption lines, which are 
fit to a single-component Gaussian curve-of-growth (COG) 
(Spitzer 1978). Although it appears from the line widths 
that the line of sight may traverse more than a single veloc- 
ity component, no high resolution optical or UV spectrum 
exists which would allow us to model a more complicated 
velocity structure. We therefore made the simplifying as- 
sumption that a single interstellar cloud with a Maxwellian 
velocity distribution is responsible for the absorption. The 
estimates of the column density (N) and Doppler param- 
eter (6) derived from the COG are likely to include sys- 
tematic errors due to this assumption. For example, the 
derived 6— value will not be a simple quadrature combi- 
nation of thermal and turbulent velocity components, as 
it would be for a true single cloud. Instead, it will also 
reflect the spread in velocities of the multiple clouds that 
are likely to lie along a sight line of this length. Jenkins 
(1986) has discussed the systematic errors associated with 
the COG technique when there are multiple clouds of var- 
ious strengths along a sight line, and has shown that if 
many non-saturated lines are used, and if the distribution 
of N and b— values in the clouds is not strongly bimodal, 
then the column density is likely to be underestimated by 

<; 15%. 

Despite these shortcomings, the COG technique has the 
virtue of being unaffected in principle by convolution with 
the instrumental line spread function (LSF). In practice, 
it must be used with caution because of the problem of 
blending with neighboring lines and, especially for a target 
like LSE 44, because it is difficult to establish the proper 
continuum in some spectral regions, presumably due to the 
presence of many unidentified stellar absorption features. 

The second technique is profile fitting (PF). For all 
species except H I we used the code Owens . f , which has 
been applied to FUSE spectra for many previous investi- 
gations. This code models the observed absorption lines 
with Voigt profiles using a \ 2 minimization procedure with 
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many free parameters, including the line spread function, 
flux zero point, gas temperature, and turbulent velocity 
within the cloud. We split each spectrum into a series 
of small sub-spectra centered on absorption lines, and fit 
them all simultaneously Each fit typically includes about 
60 spectral windows, and approximately 100 transitions of 
H I, D I, O I, N I, and H 2 ( J = to 5). The wavelengths 
and oscillator strengths of the lines used in the profile fit- 
ting analysis are given in Table 4. Note that there are fewer 
than 100 entries in the table because some lines appear in 
multiple detector segments. D I, O I, and N I were as- 
sumed to be in one component, H 2 in a second, and H I in 
a third. Velocities between windows were allowed to shift 
to accommodate inaccuracies in the FUSE wavelength cal- 
ibration. The line spread function (LSF) is fixed within a 
spectral window, but is allowed to change in a given de- 
tector segment from one window to another, and between 
segments. This is consistent with the performance of the 
FUSE spectrographs (Sahnow et al. 2000). No system- 
atic relationship between the LSF and the derived column 
density or b— value of any species was observed. 

Only one interstellar component for a given species was 
assumed along the sight line although this is unlikely to be 
true for such a distant target. Hebrard et al. (2002) dis- 
cuss tests they performed on an extensive list of potential 
systematic errors using profile fitting, including the sin- 
gle component assumption. In these tests they found that 
fits obtained with up to five interstellar components gave 
the same total column density within the 1 a error bars, 
as long as saturated lines are excluded from the analysis, 
which we have been careful to do. Thus, we report total 
integrated column densities along the sight line. 

It is possible that the temperature and turbulent ve- 
locities differ from one cloud to another. To test this we 
did profile fits with D I, O I, and N I in three indepen- 
dent components. The column density estimates differ by 
less than 1 a from those obtained by assuming a single 
component. The b— values also agree with that obtained 
by assuming a single component, although they are not 
well-constrained because only unsaturated lines are used. 
These species could also have different ionization states if 
they reside in clouds with different temperatures. How- 
ever, the ionization potentials of D I and O I and, to a 
slightly lesser extent, N I are so similar that this is un- 
likely to be a problem. 

The Fe II absorption lines detected in the spectrum are 
abnormally broad for interstellar lines, about 10 km s" 1 , 
with a temperature of ~ 2 x 10 5 K, compared to ~ 6 km s _1 
for N I and O I (see §4.2). The cause of this broadening 
is not known, but Fe II is not present in the stellar atmo- 
sphere because the temperature is too high. In the fits we 
did not include Fe II in the component containing D I, O I, 
and N I, but this has no effect on the determination of the 
interstellar column densities of these species. 

All laboratory wavelengths and oscillator strengths used 
in this work are from Morton (2003) and Abgrall et 
al. (1993a, 1993b). Neither the COG nor profile fitting 
techniques include oscillator strength uncertainties in the 
error estimates. The application of Owens . f to FUSE data 
is described in greater detail by Hebrard et al. (2002). 

4.1. The D I Analysis 



We determined N(D I) using only the profile fitting 
technique. We attempted to construct a COG, but we 
were unable to determine unambiguously the proper place- 
ment of the continuum in the presence of the strong H I 
lines separated from the D I lines by +81 km s _1 . The 
situation is particularly difficult for this sight line due to 
the unusual shape of the neighboring continuum around 
the A916.180 line, and the unaccounted for absorption on 
the red wing of the H I A920.712 line. Our tests indicated 
that failure to recognize such problems caused COG es- 
timates of N(D I) to be approximately 0.3 dex too low. 
Once the continuum is defined in the COG analysis, it 
is fixed. By contrast, although continuum placement er- 
rors will also lead to column density errors with profile 
fitting, this technique allows the shape of the continuum 
to change as part of the fitting process. The continuum 
was modeled with polynomials of order — 4 to allow for 
a large range of continuum fits. We required acceptable 
fits simultaneously to all five D I lines used in the anal- 
ysis (see below). In addition, removing any of the lines 
from the analysis did not change the final column density 
estimate significantly. This was not true for the COG fit, 
which was highly sensitive to the continuum placement at 
the A916.180 line. Removal of this line changed the COG 
result by more than 0.2 dex. Therefore, by virtue of its si- 
multaneous and self-consistent fit over multiple lines, and 
its sensitivity to the exact way in which D I blends with 
the blue wing of the neighboring H I line, profile fitting 
is the best method for measuring the deuterium column 
density. 

Five D I absorption profiles over three separate spectral 
lines were used in the profile fits: A916.180 in the SiClB 
channel, and A919.100 and A920.712 in both the SiClB 
and SiC2A channels. Longer wavelength D I lines are too 
saturated to provide additional constraints on the column 
density estimate, and the convergence of the H I Lyman 
series precludes using shorter wavelength lines. The D I 
A916.931 and A917.880 lines are too strongly blended with 
O I and H 2 lines to be used. Although the same transi- 
tion observed in separate FUSE channels potentially are 
subject to systematic errors arising from similar effects, 
such as continuum placement and blending, they do rep- 
resent independent measurements with distinct spectro- 
graphs. Thus, we include all five profiles in our analysis. 

To properly measure the interstellar D I, the effects of 
the stellar H I Lyman series and He II absorption must be 
removed. This was done in two ways. First, the stellar 
model described in §3 was shifted in velocity space and 
scaled in flux until the stellar H I damping wings matched 
the observed spectrum away from the cores of the interstel- 
lar H I lines. The resulting average velocity difference was 
AV = V* - Vism = -39 ± 5 km s" 1 and -32 ± 13 km s" 1 
for the SiClB and SiC2A spectra, respectively. This com- 
pares favorably to AV = —33 ± 6 km s _1 based on the 
measured velocities in the high-dispersion IUE spectrum 
of photospheric lines (N IV, N V, O VI, Si VI, S V, C III, 
and C IV) and interstellar lines (Si II, Si III, and S II). The 
spectral resolution of IUE in this mode is approximately 
25 km s- 1 (FWHM). 

The stellar model is only an approximation of the true 
absorption from the atmosphere of LSE 44. Since not all 
species present in the atmosphere are properly accounted 
for in the model, and since some of the atomic constants 
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included in the model are not well-known, using the model 
may in fact introduce systematic errors. To check for this 
we also applied the profile fitting without the model, us- 
ing instead up to a fourth-order polynomial to fit the stel- 
lar spectrum and continuum. This method has been used 
many times before in similar analyses; see, for example, 
Hebrard & Moos (2003). As described below, the two 
methods give column densities which agree within the er- 
rors. 

Several parameters are free to vary through the fits, in- 
cluding the column densities, the radial velocities of the 
interstellar clouds, and the shapes of the stellar continua. 
Some instrumental parameters are also free to vary, includ- 
ing the widths of the Gaussian line spread function (LSF), 
which is convolved with the Voigt profiles within each spec- 
tral window. The averaged width of the LSF that we found 
is 5.8 pixels, with a 0.8-pixel la dispersion (full widths at 
half maximum). These FUSE pixels are associated with 
CALFUSE version 3, and are 0.013A in size. Note that 
CALFUSE versions 1 and 2, which was used in most pre- 
viously published FUSE studies, produced pixels approxi- 
mately half this size. This explains why the LSF reported 
previously was approximately about 11 pixels wide (see, 
e.g., Friedman et al. 2002; Hebrard et al. 2002). Some 
examples of fits are plotted in Figure 03 

Using the stellar model we find log N(D l) sm — 15.84 ± 
0.08. Using a polynomial fit to the stellar plus interstel- 
lar continuum, we find N(D I) po i y = 15.89 ± 0.08. These 
errors include our estimates of the systematic errors asso- 
ciated with the stellar normalization. Our final estimate is 
the mean of these two values 5 , log 7V(D I)= 15.87 ± 0.08, 
where the error of the mean has not been reduced below 
the individual values since JV(D I) sm and N(D l) po iy are 
not independent estimates of the column density. 

4.2. The Ol Analysis 

The continuum placement problem is much less severe 
for the O I lines because the spectral regions adjacent to 
the lines used in the analysis are generally smooth and 
well-behaved, without the presence of very strong adjacent 
lines. Thus, we use both the profile fitting and curve-of- 
growth techniques to estimate N(0 I). The profile fitting 
analysis used only the O I A974.070 line, because this line 
is almost completely unsaturated. We did verify that the 
fit to other O I lines in the spectrum was acceptable, but 
they provide almost no additional constraint on the col- 
umn density because they are saturated. One might be 
concerned that an erroneous /—value for this transition 
would cause an inaccurate estimate of N(0 I). However, 
Hebrard et al. (2005) discuss this issue, and show that it 
is consistent with column densities derived from stronger 
O I lines in the case of BD+28°4211, and derived from 
the weak A1356 intersystem line in the case of HD 195965. 
The O I A974.070 line is blended with two H 2 lines, Ly- 
man 11-0R(2) A974.156 and Werner 2-0Q(5) A974.287, and 
absorption from these species was accounted for when cal- 
culating N(0 I), as described below. There is also an H 2 
J = 6 line at nearly the same wavelength, but the column 
density in this rotational level is too small to be of concern. 
The result of our analysis is N(0 l) pf = 17.53±g;?|. 

To construct the curve-of-growth the equivalent widths 

5 Unless otherwise noted, errors on all quantities in this paper are 2 a. 



of the O I lines were measured in the following way. First, 
low-order Legendre polynomials were fitted to the local 
continua around each O I line. The lines were integrated 
over velocity limits chosen to exclude absorption from ad- 
jacent lines. Only isolated lines were selected so that no 
deblending was required. Scattered light remaining after 
the CALFUSE 3 reduction is negligibly small, and no addi- 
tional correction was required. The measured equivalent 
widths are given in Table 5. The estimated errors from 
this method are described in detail by Sembach & Savage 
(1992). They include contributions from both statistical 
and fixed-pattern noise in the local continuum. Contin- 
uum placement is particularly difficult and subject to er- 
ror due to the presense of many unidentified stellar lines. 
To estimate the magnitude of such errors, we did trials 
in which we drew the continuum at its maximum and 
minimum plausible locations, and compared the calculated 
equivalent widths with the best estimated placement. We 
added 2 mA to the error of each measured O I line equiv- 
alent width to account for these placement errors. This 
is consistent with the magnitude of this systematic error 
determined in previous analyses of similar sight lines (e.g., 
Friedman et al., 2002). 

Due to the blending with the H 2 lines we initially did 
not include the O I A974.070 line. However, this line is 
extremely important because it is the only O I line that 
is not substantially saturated, and therefore provides the 
greatest constraint on N(0 I). We constructed curves of 
growth for the J = 2 and J — 5 rotational levels of H 2 , 
calculated column densities and b— values, determined the 
appropriate Voigt profiles, and used them to remove the 
H 2 signature from the blended line. Then the equivalent 
width of the O I line was measured in the usual way. The 
H 2 analysis is described in Section 4.5. The O I curve 
of growth is shown as the solid line in Figure H3 The 
best fit column density is log N(0 1)cog = 17.59lg;^. 
Combining this with our profile fitting result, we arrive at 

our best estimate of the O I column density, log N(0 I)= 
1 7 K 7 +o.2i 

11 - 01 -0.15- 

We note that the O I A919.917 line has been excluded 
from the analysis. We were unable to get a good fit of this 
line with profile fitting, and it fell well below the curve- 
of-growth. This is also true for Feige 110 (Hebrard et 
al. 2005), whose spectrum is similar in many respects to 
that of LSE 44, as discussed in §3. It is unlikely that this 
line has a significantly incorrect /—value since the absorp- 
tion from this line is consistent with absorption from other 
O I lines in, for example, the spectrum of BD+28°4211 
(Hebrard & Moos 2003). It is possible that there is uniden- 
tified absorption adjacent to this line, which depresses the 
local continuum, causing an underestimate of the equiva- 
lent width of this line. Alternatively, the location of this 
line on the COG may be an indication that there are multi- 
ple components along the sight line, with different column 
densities. This would cause a "kink" in the COG, which 
might be more obvious if additional O I lines were available 
to include on the plot. However, with no high resolution 
spectra available, we are unable to confirm this hypothesis. 
If we include the A919.917 line, the single-component COG 
gives log N(0 i)= I7.52±g:|| and b = 6.33±£;^, which is 
shown as the dashed line in Figure This differs from our 
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best estimate column density by less than la. 

It is worth emphasizing the importance of non-saturated 
lines in the determination of N(0 I) even at the expense of 
the potential introduction of systematic errors associated 
with using strongly blended lines. Computing a COG ex- 
cluding the O I A974.070 (and A919.917) lines gives N(0 I) 
= 17.221q'29, about 2a less than the result obtained when 
including the weak line. This is shown as the dotted line 
in Figure^] The fit to the remaining spectral lines is quite 
acceptable, but is clearly inconsistent with the weak line. 
Indeed, this line was not used in the initial analysis of 
Feige 110 (Friedman et al. 2002). We have recently in- 
cluded A974.070 in a re-analysis of iV(0 I) along this sight 
line. The COG and profile fitting methods are now in good 
agreement, and the column density estimate has increased 
by 0.33 dex compared to the original estimate. This is 
discussed more fully in Hebrard et al. (2005). 

4.3. The Nl Analysis 

As was the case with D I, only profile fitting was used to 
determine N(N I). The COG technique was not used be- 
cause of the uncertainties in establishing continuum levels 
on each side of the N I lines. 

Three lines were used in the profile fitting analysis: 
A951.079, A951.295, and A955.882. The method is the 
same as was used for the deuterium analysis, and is de- 
scribed in §4.1. The result is log N(N I) = 16.43 ± 0.14. 

4.4. The HI Analysis 

In the direction of LSE 44, the Lya line provides the 
best constraint on the total H I column density because its 
radiation damping wings are very strong. However, Lya 
is not accessible to FUSE. Ly/3 displays modest damping 
wings, and we attempted to model the absorption profile 
in order to estimate N(H I). We were unable to do this 
primarily because of the presense of two H2 lines, J=0 
A1024.230 and J=2 A1026.526, on the blue and red wings 
of the Ly/3 absorption, respectively. As discussed in §4.5, 
all lines of H2 in the J=0 rotational level lie on the flat 
portion of the COG, so we are unable to make even a 
rough estimate of the column density. The situation for 
J=2 is only slightly better, with the weakest lines still con- 
siderably saturated. Since these lines are located almost 
exactly at the regions of the Ly/3 profile which are most 
sensitive to constraining a model profile, it is impossible 
to obtain a good estimate of N(H I). The higher order 
Lyman lines in the FUSE spectrum are even weaker rel- 
ative to the stellar absorption than is Lya or Ly/3. Cor- 
recting these lines for stellar absorption would therefore 
require more precise knowledge of the proper flux scaling 
of the model, and of the velocity difference between the 
stellar and interstellar components. We have examined 
the range of b— value/column density combinations allowed 
by the higher Lyman series lines in the FUSE spectrum, 
along with the uncertainties in correcting for the stellar 
absorption, and we find that these lines do not constrain 
A(H I) with sufficient precision to provide an interesting 
D/H measurement. No suitable HST spectra of the LSE 
44 Lya line have been obtained. Consequently, we have 
used the only high-dispersion echelle observation of this 
star obtained with the short wavelength prime camera on 
IUE (exposure ID SWP29086) to estimate A(H I) toward 



LSE 44. This IUE mode provides a resolution of ~25 
km s _1 FWHM and covers the 1150-1950 A range. 

The IUE observation was obtained on 1986 August 29 
with an exposure time of 16.5 ksec. We retrieved the 
IUE spectrum from the STScI MAST archive with both 
IUESIPS and NEWSIPS processing. Since there are sev- 
eral known problems with NEWSIPS processing of high- 
dispersion observations (Massa et al. 1998) that can ad- 
versely affect the A(H I) measurement, particularly the 
background subtraction in the vicinity of Lya where the 
orders are closely spaced on the detector (Smith 1999), we 
generally prefer IUESIPS spectra. Nevertheless, we com- 
pared the IUESIPS and NEWSIPS versions of the spec- 
trum at the outset. In addition to the standard processing, 
we applied the Bianchi & Bohlin (1984) correction method 
to the IUESIPS spectrum to compensate for scattered light 
from adjacent orders in the background regions near the 
interstellar Lya line. 

4.4.1. Analysis Method 

To measure the total interstellar H I column density, 
we used the method of Jenkins (1971); see also Jenk- 
ins et al. (1999) and Sonneborn et al. (2000). In 
brief, we constrained the H I column density using the 
Lorentzian wings of the Lya profile, which have optical 
depth r at wavelength A given by r(A) = N(R I)cr(A) = 
4.26 x 10- 20 iV(H I) (A - A )- 2 . Here A is the centroid of 
the interstellar H I absorption in this direction, which was 
determined from the N I triplet at 1200 A. We estimated 
the value of N(H I) that provides the best fit to the Lya 
profile by minimizing % 2 using Powell's method with five 
free parameters: (1) N(H I), (2)-(4) three coefficients that 
fit a second-order polynomial to the continuum (with a 
model stellar Lya line superimposed, see below), and (5) 
a correction for the flux zero level. We then increased (or 
decreased) N(H I) while allowing the other free param- 
eters to vary in order to set upper and lower confidence 
limits based on the ensuing changes in \ 2 ■ 

In addition to the Lya absorption profile due to inter- 
stellar H I , a subdwarf O star like LSE 44 will have a sub- 
stantial stellar Lya absorption line with broad, Lorentzian 
wings as well. Neglect of this stellar line when setting the 
continuum for fitting the interstellar Lya will lead to a 
substantial systematic overestimate of A(H I). We have 
used the stellar atmosphere models discussed in §3 to ac- 
count for the stellar Lya line. The offset of the stellar line 
centroid with respect to the interstellar lines was deter- 
mined by comparing the velocities of several stellar and 
interstellar lines in the FUSE spectrum, and the uncer- 
tainty in the determined offset has a negligible impact on 
the derived N(R I). We use the most shallow and deepest 
stellar profiles derived from a grid of models, as described 
in §3, to set upper and lower confidence limits on 7V(H I) 
in the following section. 

4.4.2. N(Hl) Results 

Figure[7]shows the IUE spectrum of LSE 44; both panels 
plot the same data but the upper panel shows a broader 
wavelength range to enable the reader to inspect the fit 
to the continuum away from the strong Lya absorption. 
Overplotted on the spectrum with dotted lines are the fits 
that provide 2 a upper and lower limits on the interstellar 
iV(H I). The deepest model stellar Lya line is assumed for 
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the lower limit on N(H I), and the most shallow stellar 
Lya model is assumed for the upper limit. From these fits 
we derive log 7V(H I) = 20.52+g;|g (2 a). 

There are, of course, many potential sources of system- 
atic error in the N(H I) measurement. We have already 
discussed the stellar Lya line. Strong, but nevertheless 
unrecognized, narrow stellar lines within the Lya profile 
are another potential source of systematic error. We have 
included stellar lines in the continuum regions used to es- 
timate a for the \ 2 calculation, but nevertheless future 
higher resolution and signal-to-noise observations that re- 
solve strong stellar lines might yield a substantially dif- 
ferent result. However, it seems likely that the greatest 
source of systematic errors in this case is location of the 
zero-point flux. In a saturated line the core should be 
broad and completely black. However, examination of Fig- 
ure [7| shows that the core has significant undulations. We 
have decided that the flattening of the spectrum between 
approximately 1216. 5A and 1217. 8A gives the best indica- 
tion of the zero point. The spectral regions used in the fit 
to determine N(H I) are indicated in Figure[7]by the dark 
horizontal lines below the spectrum. 

We recognize that there are other plausible zero-point 
placements. For example, Friedman et al. (2005) initially 
used the spectral region indicated by the gray horizontal 
line in Figure [7] to determine the zero point, completely 
excluding the flat part in the red side of the airglow emis- 
sion. The result of this fit was log 7V(H I) = 20.41±°;^. 
However, we ultimately decided that the flattening to the 
red of the airglow line better represents the expected shape 
of the line core, and the undulations on the blue side make 
this region an unreliable indicator of the zero point. The 
undulations could have many sources, including improper 
joining of spectral orders in the IUE echellogram. We note 
that if the region on the blue side is used to set the zero 
point, we obtain a lower value for N(H I), which in turn 
makes D/H even higher than the value we obtain in the 
next section. 

4.5. The H2 Analysis 

Accurate measurements of the column densities of 
molecular hydrogen in various rotational levels are impor- 
tant so that the H2 lines may be deblended in order to 
measure the column densities of D I and metals. For ex- 
ample, as noted in §4.2, the weak O I A974.070 line is 
blended with J = 2 and J = 5 lines. We attempted to 
measure N(R2) for J = — 5 so that we could estimate 
the total H2 column density and determine the molecu- 
lar fraction of the gas along the sight line. However, the 
J = and J = 1 rotational levels could not be accurately 
determined. In each case there are no lines strong enough 
to show well-developed damping wings, and no lines weak 
enough to be on the linear portion of the curve-of-growth. 
The situation is slightly better for J = 2, which has some 
weaker lines but none that are in the true linear regime. 
We constructed a COG, but the column density is not 
well-constrained. For J — 3 — 5, weaker lines are available 
and the errors are lower. The measured equivalent widths 
are given in Table 6, and the column densities are given in 
Table 7. The H 2 COGs are shown in Figure El 

5. RESULTS AND DISCUSSION 



We have measured the column densities of D I, O I, N I, 
H I, and H2 along the sight line to LSE 44. This target is of 
particular interest because it is one of only five for which 
such values are known for Galactic targets at distances 
greater than about 500 pc and hydrogen column densities 
log iV(H l)> 20.5. Table 8 gives a summary of the results 
of this study, and Figure shows our D/H value compared 
to previous measurements. Our principle results are D/H 
= (2.24±i;||) x 10- 5 , D/O = (1.99tJ;jj?) x 10~ 2 , and O/H 

= (l.Uton) x 10~ 3 - Of the targets with the H I col- 
umn densities only LSE 44 and PG0038+199 (Williger et 
al. 2005) have high D/H, about 2.2 x 10~ 5 . D/H for the 
remaining targets are all < 1.0 x 10 -5 . 

The spatial variation of D/H has drawn increased at- 
tention in recent years. The large number of measure- 
ments made using FUSE, together with earlier results from 
IMAPS, HST, and Copernicus, led to initial speculation 
that the observed variations in D/H and O/H were due 
to differing levels of astration (Moos et al. 2002). Mod- 
els of galactic chemical evolution (dc Avillez & Mac Low 
2002; Chiappini et al. 2002) attempt to account for the 
spatial variability in terms of mixing scale lengths and 
times. One of the greatest difficulties for such models 
is to properly account for the astration factor required 
to reduce the primordial deuterium abundance to that 
observed along some sight lines in the Galaxy. The pri- 
mordial value, estimated from gas-phase measurements of 
low-metallicity material at high redshift toward 5 quasars, 
is (D/H) p „ m = 2.78to.H x 10 ~ 5 (Kirkman et al. 2003). 
This is consistent with the value inferred from the ampli- 
tude of the acoustic peaks in cosmic background radiation 
(CBR) measurements made with WMAP, (D/R)cbr = 
2 - 37 ±o:2i x 10 ~ 5 (Spergel et al. 2003). Thus, astration 
factors of approximately 1.6 or 1.9 are required to account 
for Local Bubble D/H abundances determined by Wood 
et al. [(D/H)iB = (1-56 ± 0.04) x 10~ 5 ], and Hebrard 
& Moos [(D/H) LB = (1.32 ±0.08) x 10~ 5 ], respectively. 
An astration factor of approximately 3 is required to ac- 
count for the mean D/H value of the four highest column 
density sight lines in Figure Hfb), and this is pushing the 
limits of most chemical evolution models (Tosi et al. 1998). 
These models will also have to account for the deuterium 
abundance measured in the Complex C, a low-metallicity, 
high-velocity cloud falling in to the Milky Way. Sembach 
et al. (2004) found (D/R) complexC = (2.2 ± 0.7) x lO" 5 
and (D/O) coropte c = 0.28 ± 0.12 (Lr). 

Hebrard & Moos (2003) suggest that the low D/H value 
at large distances truly reflects the present-epoch D/H 
abundance. This conclusion is primarily based on D/O 
and O/H measurements. If this is correct, it presents a 
challenge to models of galactic chemical evolution, which 
not only have to account for high astration values, but also 
the observed spatial variation of these quantities. 

Jura (1982) first proposed that deuterium might be de- 
pleted onto dust grains in the interstellar medium. Draine 
(2004a, 2004b) extended this by suggesting that D/H in 
carbonaceous dust grains can exceed the overall D/H ratio 
by a factor greater than 10 4 , if the material is in thermody- 
namic equilibrium. Thus, it is possible that a significant 
reservoir of deuterium is depleted onto these grains, as 
long as the metallicity is approximately solar or greater. 
Prochaska, Tripp, & Howk (2005) have tested this idea by 
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examining the correlation between D/H and the highly de- 
pleted ion Ti II. They find these are correlated at greater 
than 95% confidence level. They emphasize, however, that 
since not all data points support this relationship, mea- 
surements of distant targets, such as the one in the current 
study, are most needed to further constrain this relation- 
ship. At present, we are unable to comment directly on 
this hypothesis because the Ti II abundance along this 
sight line is not known. 

Wood et al. (2004) and Linsky et al. (2005, in prepa- 
ration) have applied Draine's suggestion to the obser- 
vational data, and suggest the presence of three D/H 
regimes, as shown in Figure |§1 In the Local Bubble, 
where log 7V(H I)<^ 19.2, most of the deuterium is in the 
gas phase, presumably because of the recent passage of 
strong shocks from supernovae or stellar winds. Thus, 
(D/R) LBg = (1.56 ± 0.04) x 10~ 5 is relatively high, spa- 
tially uniform, and approximates the true local Galactic 
disk value. At great distances, where log N(R I)> 20.5, 
sight lines traverse many regions, some recently shocked 
and others not. D/H is likely to be indicative of an aver- 
age of these regions. However, the highest column density 
sight lines, which tend to be the lengthiest ones, will be 
weighted toward the densest, coldest regions, and presum- 
ably the most heavily deuterium-depleted clouds. This 
may explain the low D/H values toward the majority of 
the most distant targets. In the intermediate column den- 
sity regime, 19.2 <; log(A^(HI)) <; 20.5, sight lines tra- 
verse fewer regions with different shock event histories and 
therefore exhibit more variability. 

Hebrard & Moos (2003) show that D/O tends to be 
uniform within the Local Bubble, with a mean of (3.84 ± 
0.16) x 10~ 2 (1<t), and variable beyond it. D/N is similar, 
but exhibits somewhat more variability due to ionization 
effects. At the highest column densities D/O and D/N 
converge to mean values that are lower than the Local 
Bubble values. Indeed, for all seven sight lines for which 
log N(U I) > 20.1 and log N(D I) > 15.4 (Hoopes, et al. 
2003; Wood et al. 2004; Williger et al.), including LSE 44, 
D/O lies in the relatively narrow range (1.8 — 2.6) x 10~ 2 . 
The fact that D/O is much more uniform than D/H is a 
surprise, and suggests to some investigators that measure- 
ments of N(H I) may suffer from additional unidentified 
errors (Hebrard et al. 2005). 

The ~ 100 pc length scale in the Local Bubble over 
which D/H is roughly constant (Figure EJ) may be a nat- 
ural consequence of supernova-driven mixing within the 
ISM. Moos et al. (2002) discuss a simple model which has 
two time scales: the mixing time and the time between 
supernovae (SNe). In a fully ionized ISM at ~ 10 6 K, the 
adiabatic sound velocity is c s ~ 150 km s -1 . For a spher- 
ical parcel of gas with radius 100 pc, the crossing time is 
t s w 7 x 10 5 years. Other mechanisms, such as the Alfven 
crossing time for a cloud with a magnetic field, are gener- 
ally slower. 

A supernova rate of 0.02 yr _1 (Cappellaro et al. 1993) 
spread uniformly over the galactic disk of radius 8.5 kpc, 
yields approximately 8.8 x 10~ n SNe pc~ 2 yr _1 . In the 
100 pc radius parcel of gas considered above, there would 
be one supernova every 4 x 10 5 yr. Therefore, in this sim- 
ple model, tsN ~ is- Thus, one might expect that within 
the Local Bubble, the material is well- mixed, and D/H, 
O/H, and N/H will be approximately constant, with rel- 



atively small target-to-target variations. Toward targets 
at greater distances, different supernovae events may have 
given rise to interstellar regions with different enrichment 
levels, and these do not have time to mix thoroughly with 
each other. See Moos et al. (2002) for additional dis- 
cussion of this model. More sophisticated treatments of 
mixing, including the effects of supernovae and infall of 
material into the Galactic disk, are given by de Avillez 
(2000), de Avillez & Mac Low (2002), and Chiappini et al. 
(2002). 

If this model is correct, it is hard to explain the uni- 
formity of D/O for the high column density sight lines, 
unless a different mechanism forces this to be the true 
Galactic value. In addition, for these targets D/H and 
O/H vary by factors of 2.9 and 4.6, respectively, but there 
is no convincing anticorrelation between them, as might 
be expected if O is produced as D is destroyed in stellar 
interiors. Steigman (2003) reported a tentative detection 
of such an anticorrelation based on IMAPS and the ini- 
tial set of FUSE results, all of which had column densities 
lower than those considered here. 

Neither the depletion hypothesis nor the low Galactic 
D/H hypothesis appears to be consistent with all of the 
data currently available. Hebrard et al. (2005) show 
that the dispersion of D/H is greater than that of D/O 
in the sense that there are no high D/O values at large 
distances or large column densities. They argue against 
the depletion hypothesis since for sight lines beyond the 
Local Bubble, the uniformity of D/O and O/H (Meyer, 
Jura, & Cardelli 1998; Meyer 2001; Andre et al. 2003) 
strongly imply that D/H must be uniform as well, unless 
D/O and O/H are correlated in exactly the same way. 
They emphasize that measurements of D/H may be sub- 
ject to fewer systematic errors then measurments of D/O 
and O/H. On the other hand, nobody has identified an 
error that could account for the D/H variability shown in 
Fig-El Two observational tests will help clarify these inter- 
pretations. First, additional measurements of the deple- 
tion of refractory species, such as Fe, Si, or Ti (Prochaska, 
Tripp, & Howk 2005) along these sight lines would re- 
veal whether a correlation exists between depletion and 
the gas phase abundance of D/H. Such a correlation would 
provide strong evidence for the depletion hypothesis. Sec- 
ond, Lecour et al. (2005) measured N(HD) and N(H 2 ) 
along highly reddened sight lines, and inferred (D/H) = 
(2.0 ± 1.1) x 10~ 6 an extremely low value. Linsky et al. 
(2005, in preparation) suggest that this is due to high lev- 
els of deuterium depletion, and they discuss methods to 
investigate this which would test the depletion hypothe- 
sis. 

It is possible that the systematic errors in these mea- 
surements are larger than the investigators have estimated, 
leading to erroneous conclusions about the variability of D, 
O, and N. This may also be the cause of the apparent dis- 
persion in D/H in high redshift environments (Crighton 
et al. 2004), where very little variability is expected since 
the gas has experienced little astration and has low metal- 
licity. A potential alternative explanation of the variabil- 
ity would be a mechanism of non-primordial deuterium 
production (Epstein, Lattimer, & Schramm 1976). Many 
mechanisms have been examined carefully including, for 
example, D production in stellar flares (Mullan & Linsky 
1999), but none has been shown to produce enough to ma- 
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terially change the composition of the ISM (Prodanovic & 
Fields 2003). See Lemoine et al. (1999) for a review of a 
variety of potential production mechanisms. 

Combining results from Hebrard & Moos (2003), Wood 
etal. (2004), and Hoopeset al. (2003), it is seen that D/H, 
O/H, and N/H are factors of 1.5 — 3 higher toward LSE 
44 than toward most other stars with log JV(H I) > 20.5. 
This at least suggests the possibility that our estimate 
of A^(H I) may be underestimated by approximately this 
amount. We discussed in §4.4 the problems of determining 
the proper flux zero point in the core of the H I line. If 
our zero-point determination is correct, then our error in 
N(H I) is just consistent with these abundances being a 
factor of 1.5 greater than our best estimate, but not with 
a factor of 3. Unfortunately, without improved measure- 
ments of the Lya profile, better estimates of iV(H I) are 
not likely to be obtained. 

Studies by the FUSE deuterium team have been con- 
centrating on distant targets in order to probe more com- 
pletely both the transition and high column density re- 
gions displayed in Figure ED LSE 44 is one of several 



new targets in this regime. Additional observations of tar- 
gets at these large distances will be required to determine 
whether depletion has caused the observed spatial varia- 
tions in the gas phase D/H ratio, or if the Galactic disk 
value is really below 10 -5 and different histories of as- 
tration are of primary importance, or if some unknown 
non-primordial source of deuterium is responsible for the 
variability. 

We thank Dave Kilkenny for kindly providing the opti- 
cal spectrum of LSE 44, Gerry Williger for reducing the 
optical data, and Brian Wood for providing the data in 
Figure |§1 in numerical form. This work is based on data 
obtained for the Guaranteed Time Team by the NASA- 
CNES-CSA FUSE mission operated by The Johns Hopkins 
University. This work used the profile fitting procedure 
Owens . f developed by M. Lemoine and the French FUSE 
team. Financial support has been provided by NASA 
contract NAS5-32985. SDF is supported by NASA grant 
NNG04GH19G and TMT by NASA grant NNG04GG73G. 
GH is also supported by CNES. 
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Table 1 
LSE 44 Observation Log 



Instrument 


Date 


Program ID 


Exposures' 1 


Exposure Time (ks) 


FUSE 


24 April 2002 


P2051602 


14 


13.9 


FUSE 


09 August 2003 


P2051603 


9 


6.2 


FUSE 


10 August 2003 


P2051604 


14 


5.6 


FUSE 


11 August 2003 


P2051605 


25 


18.2 


FUSE 


12 August 2003 


P2051606 


12 


8.5 


FUSE 


03 February 2004 


P2051607 


20 


15.7 


FUSE 


04 February 2004 


P2051608 


14 


11.1 


FUSE 


05 February 2004 


P2051609 


11 


6.3 


IUE 


29 August 1986 




SWP 29086 


16.5 



a Not all exposures for each visit were included in the analysis of the FUSE data. 



Table 2 
Exposure Time by Channel 



Channel 


Exposure Time (ks) 


LiFIA 


64.5 


LiFIB 


64.8 


L1F2A 


50.6 


L1F2B 


51.3 


SiCIA 


61.2 


SiClB 


61.2 


SiC2A 


61.9 


SiC2B 


60.4 



Table 3 

Target and Sight Line Parameters for LSE 44 



Quantity 


Value a 


Reference 13 


Spectral Type 


sdO 


1 


(l,b) 


(313?37,+13?49) 


1 


d (pc) 


554 ± 66 


2 


z (pc) 


129 ±15 


2 


V 


12.45 


1 


U -B 


-1.18 


1 


B-V 


-0.24 


1 


Eb-v 


0.05 ±0.03 


2 


T cff (K) 


38700 ± 1000 


2 


\ogg (cm s 2 ) 


5.5 ±0.1 


2 


log [N(Re)/N(R)] 


-2.8 ±0.1 


2 



a Errors are la. 

References: (1) Drilling 1983; (2) this study. 
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Table 4 

Spectral Lines Used in the Profile Fitting Analysis 



G 

opGCIGS 




J V ell UL^ 




Wavplpno-fVi ( K \ 
V V d V cltllti LI1 In. 1 


/ V ell UL^ 


HI 


919.351 


1.20E-03 


H 2 J=3 


1028.985 


1.74E-02 


HI 


920.963 


1.61E-03 




1043.502 


1.08E-02 


D I 


916.18 


5.78E-04 




1056.472 


9.56E-03 


D I 


919.101 


1.20E-03 




928.436 


3.30E-03 


D I 


920.712 


1.61E-03 




933.578 


1.95E-02 


NI 


951.079 


1.69E-04 




934.789 


7.09E-03 


NI 


951.295 


2.32E-05 




935.573 


6.57E-03 


NI 


955.882 


5.88E-05 




967.673 


2.28E-02 


01 


974.07 


1.56E-05 




970.56 


9.75E-03 


Fe II 


1055.262 


6.15E-03 




985.962 


8.28E-03 


Fe II 


1062.152 


2.91E-03 


H 2 J=4 


1023.434 


1.05E-02 


Fe II 


1125.448 


1.56E-02 




1032.349 


1.71E-02 


Fe II 


926.897 


5.66E-03 




1044.542 


1.55E-02 


Fe II 


935.518 


2.56E-02 




1057.38 


1.29E-02 


H 2 J=0 


938.465 


9.51E-03 




919.048 


9.30E-03 


H 2 J=l 


955.708 


4.23E-03 




920.834 


1.31E-02 


H 2 J=2 


1026.526 


1.80E-02 




933.788 


1.08E-02 




1053.284 


9.02E-03 




935.958 


1.95E-02 




920.241 


1.68E-03 




938.726 


6.38E-03 




927.017 


2.33E-03 




940.384 


1.99E-03 




932.604 


4.78E-03 




955.851 


1.69E-03 




940.623 


6.03E-03 




962.151 


8.75E-03 




941.596 


3.37E-03 




968.664 


1.26E-02 




974.156 


1.32E-02 




970.835 


1.60E-02 




975.344 


6.64E-03 




971.387 


3.51E-02 




984.862 


8.17E-03 


H 2 J=5 


994.227 

1052.496 

1061.697 

1065.596 

916.101 

928.76 

955.681 

974.286 


1.38E-02 
1.11E-02 
1.26E-02 
9.90E-03 
2.38E-03 
3.91E-03 
2.76E-02 
3.52E-02 



Table 5 
O I Equivalent Widths 3 - 



W x (ml) 



A(A) 


log A/ 


SiClB 


SiC2A 


LiFIA LiF2B 


925.446 


-0.484 


74.4 ±6.0 


71.0 ±4.5 




971.738 


1.128 


112.3 ±6.3 


105.2 ±6.4 




974.070 


-1.817 


33.3 ±4.8 


34.3 ±4.5 




1039.230 


0.974 






122.0 ±3.1 120.7 ±3.8 



a Oscillator strengths from Morton 2003. Equivalent width errors are la. 



Table 6 
H 2 Equivalent Widths 3, 





A(A) 




W x (: 


ml) 


Rotational Level 


log A/ 


SiClB 


S1C2A 


J=2 


920.241 


0.190 


57.7 ± 3.9 


53.6 ±5.4 




927.017 


0.334 


60.7 ±4.0 


61. 4± 2.2 




932.604 


0.649 


69.4 ± 2.1 


74.3 ± 2.3 




940.623 


0.754 


74.0 ± 2.6 


77.9 ±5.5 




975.344 


0.812 


75.4 ± 2.5 


77.0 ± 1.7 


J=3 


933.578 


1.260 


82.8 ± 2.3 


77.5 ± 2.1 




934.789 


0.821 


65.9 ± 2.2 


61.0 ± 2.4 




936.854 


0.283 


36.8 ± 2.4 


38.9 ±3.9 




951.672 


1.081 


77.2 ± 2.3 


77.5 ± 1.9 




960.449 


0.676 


56.1 ± 2.3 


56.1 ± 2.3 


J=4 


933.788 


1.002 


11.3 ±2.7 


11.0 ±2.1 




968.664 


1.086 


13.4 ±2.1 


14.7 ±2.6 




970.835 


1.192 


19.2 ±2.2 


19.2 ±2.4 


J=5 


938.909 


1.265 


14.7 ±4.6 


16.7 ±5.7 




940.882 


0.927 


12.8 ± 1.7 


11.5 ± 1.7 




942.685 


0.769 


7.1 ± 1.8 


6.9 ±2.5 




958.009 


0.964 


14.2 ±1.9 






974.884 


1.142 


11.4 ±2.8 


11.9 ±2.7 




983.897 


1.094 


14.8 ± 1.6 





a Equivalent width errors are la. 



Table 7 

Molecular Hydrogen Column Densities 



Rotational Level 


log N (cm" 2 ) 3 


J = 2 


- 16.20 


J = 3 


15.52 ±0.13 


J = 4 


14.27 ±0.08 


J = 5 


14.19 ±0.11 



a Errors are 2a. 



Table 8 

Summary of Results of This Study 



Quantity Value (2a errors) 



log N(D I) 




15.87 ±0.08 




log N(0 I) 




17 C7+0.21 

-0.15 




log N(N I) 




16.43 ±0.14 




log N(R I) 




20 ^0+0-20 

ZU - OZ -0.36 




D/H 


(2. 


■ 24 ±i'll) x 10" 


-5 


O/H 


(1 


.13±g-f?) x 10" 


-3 


N/H 


(8 


13^ °9) x 10- 


-5 


D/O 


(1 


•99±oi?) x 10- 


-2 


D/N 


(2 


■75±J:1|) x 10- 


-1 



13 




Fig. 1. — The SiClB spectrum of LSE 44. The positions of interstellar D I, O I, N I, and H2 absorption lines are indicated. Many of the 
unmarked lines are due to metals arising in the atmosphere of this subdwarf O star. The data have been binned by 2 pixels (~ 0.026 A) for 
display purposes only in this figure. Note that the flux scale is not the same on all panels. 
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Fig. 2. — Black line: normalized optical spectrum of the sdO star LSE 44. Red line: Best model matching the optical spectrum 
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Fig. 3. — Position of the sdO star LSE 44 on the T c g-g plane. Open circles arc parameters of observed sdB and sdOB stars. The position 
of the sdOB star Feige 110 is also indicated (Friedman et al. 2002). Solid lines are post-EHB evolutionary tracks from B. Dorman (1999, 
private communication). The numbers indicate the mass of a star in solar masses that corresponds to an evolutionary track. The hashed 
region is the ZAEHB. 
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- LSE 44 

°' 2 ~ Feige 110 

Q L I I I I I I I I I I I I I I I I I I I 

1154 1156 1158 1160 1162 

Wavelength (A) 

Fig. 4. — Portions of the FUSE spectrum of LSE 44 (red line) compared to the FUSE spectrum of the sdOB star Feige 110 (black line). 
Many unidentified lines appear in spectra of both stars. 
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S ;a l 




973.6 973,8 9/4.G 9/4.2 9/4.4 1044.0 1044.2 1044.4 1044.6 1044.8 



waselengt (A) 

Fig. 5. — Examples of spectral windows used in the profile fitting analysis. Histogram lines are the data, heavy black lines are the fits, 
light black lines are the continua, and the dashed lines are the fits for the individual species. The dotted lines are the model profiles prior to 
convolution with the LSF. The H2 lines of levels J = 2 and J = 5 are designated as H22 and H25. The X-axis is wavelength in A and the 
Y-axis is flux in units of ergs cm~ 2 s" 1 A" 1 . Eight spectral windows are displayed here, including 18 transitions. However, the complete fit 
includes typically 60 spectral windows in all FUSE detector segments, and approximately 100 transitions. 
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-3.5 



-4.0 



m 
o 



-4.5 - 



-5.0 




10 



log (fA) 



Fig. 6. — The single component, Gaussian curve of growth for O I. The best fit solution, which excludes the A919 line, is log 7V(0 I) 



17.59 



+0.13 



-0.14 



and b = 6.33 



+0.30 
0.26 



(solid line). Including all lines gives log N(0 I) 

17.22+g;|| and 6 = 



= 17.52] 
6.85 H 



and 



'•33~^o 4 4 (dashed line). The solution if 
which demonstrates that the lack of 



,+0.22 
-0.23 

both the A974 and A919 lines are excluded is log N(0 I) = 17.22^0 29 and b = 6.85± ;!^ (dotted line), 
optically thin lines can result in a large error in the column density estimate. See §4.2 for additional discussion. For clarity, at each wavelength 
the data points derived from each of the FUSE channels have been slightly separated. The inset shows the log(7V)/6— value error contours for 
the best fit solution. 
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-0.4 hiiiiiiiiiiiiiiiiiiiiiiiiiiiii - 

1200 1205 1210 1215 1220 1225 1230 

Observed Wavelength (A) 

Fig. 7. — High-dispcrsion IUE observation of LSE 44 (histogram) with the best fit H I profile (black line through spectrum) and upper and 
lower bounds (dotted lines) on N(H I) at the 2 a confidence level. The highest curve is the stellar model, described in §3. We have used the 
spectral region indicated by the dark horizontal line in the fit. Note that this uses the region from approximately 1216.5 — 1218. lA, in the 
core of the saturated line, to establish the flux zero point. We initially used the spectral region indicated by the gray horizontal line for the 
fit, but determined that this incorrectly set the zero point of the flux. See §4.4.2 for additional discussion. The upper and lower panels plot 
the same data and fits on different scales to enable the reader to inspect the details in the Lya profile as well as the continuum fit well away 
from the interstellar H I line. The emission line in the center of the saturated Lya core is geocoronal Lya emission, and is excluded from the 
analysis of the interstellar H I column density. 
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7.0 7.5 8.0 8.5 9.0 7.0 7.5 8.0 8.5 9.0 

log (Nf A) 

Fig. 8. — Curves-of-growth for H2 J = 2 — 5 rotational levels. All curves are ploted over the same range of log (NfA) to clearly display the 
rotational levels for which N(H2) is well-constrained. For J = 2 it is not possible to detemine a meaningful column density error because all 
lines are saturated. The b— value is well constrained only for J = 3, for which b = 7.54^q ^ km a -1 . We set b = 7.54 when computing the 
COG for J = 4 and J = 5, although the column densities are insensitive to the choice of b. Filled symbols are data from SiClB channel, and 
open symbols from SiC2A. 
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Fig. 9. — D/H as a function of (a) distance and (b) H I column density. In this figure, adapted from Wood et al. (2004), the la errors 
are plotted. The filled square shows the data point from this study. The arrows in panel (b) indicate the Wood at al. boundary of the Local 
Bubble, at log N(H I) = 19.2, and the start of the third D/H regime, at log N(H I) = 20.5, where D/H is averaged over many regions, leading 
to a roughly constant value. The horizontal lines indicate the average values of D/H for log N(H I)< 19.2 and log N(H I)> 20.5 (excluding 
LSE 44), respectively. Mixing caused by supernovae may be responsible for the uniform D/H distribution within ~ 100 pc. Also included in 
these plots are results from Williger et al. (2005) and Hebrard et al. (2005). 



